Ensembles of plasmonic nanospheres at optical frequencies and a problem of negative 

index behavior 
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Arrays of metallic nanoparticles support individual and collective plasmonic excitations that con- 
tribute to unusual phenomena like surface enhanced Raman scattering, anomalous transparency, 
negative index, and subwavelength resolution in various metamaterials. We have examined the elec- 
tromagnetic response of dual Kron's lattice and films containing up to three monolayers of metallic 
nanospheres. It appears that open cubic Kron's lattice exhibits 'soft' electromagnetic response but 
no negative index behavior. The close-packed arrays behave similarly: there are plasmon resonances 
and very high transmission at certain wavelengths that are much larger than the separation between 
the particles, and a 'soft' magnetic response, with small but positive effective index of refraction. It 
would be interesting to check those predictions experimentally. 

PACS numbers: 78.20.Ci, 42.30.Wb, 73.20.Mf, 42.25.Bs 



I. INTRODUCTION 

Media with strong dispersion of the refractive index 
may support backward waves, as is obvious from the fol- 
lowing relation established by Lord Rayleigh in 1877: 
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where Ug is the group index for waves with wavelength 
A, phase velocity is Vg = c/ug, group velocity v = c/n, 
where c is the light speed in vacuo. Hence, the group 
velocity may become negative in a system with large pos- 
itive dispersion. The structures that support backward 
waves have been known since early 1900s and widely used 
in antenna and electronics technologies since 1950s. Me- 
dia where the directions of phase and group velocities are 
opposite are known to produce negative refraction [l|. It 
was noted in Ref. [l| that e.g. dispersion of light close 
to excitonic frequencies in solids can be negative. It is 
also true of artificial metamaterials with strong spatial 
dispersion, like photonic crystals 0, 0] . Pafomov [4] and 
Veselago @ showed that backward waves can propagate 
in isotropic medium with simultaneously negative per- 
mittivity e and permeability fi, exhibiting negative re- 
fraction, inverse Doppler and Vavilov-Cherenkov effects. 
Currently, such metamaterials are called negative index 
metamaterials (NIM). Veselago has noticed that a slab 
of NIM with thickness d would work as a flat lens. The 
lens could not produce an image of a distant object, since 
it cannot focus a parallel beam of light, but produces a 
replica of the object if it is placed less than distance d 
away from the nearest surface of the slab. 

Recently, unusual properties of NIM systems attracted 
considerable attention followed first theoretical 0| 
and then experimental Q demonstration of feasibility 
of negative refraction by periodic metallic metamaterials 
like periodic wire meshes producing negative e [6| and 
split-ring resonators (SRR) giving negative /.t 0] in mi- 
crowave region of incident radiation. Pendry has showed 
theoretically that ideal Veselago lens can produce sub- 



wavelength resolution (perfect lens) [9l]. The effect ap- 
pears because ideal NIM slab supports surface plasmon 
modes that are in resonance with incident radiation at 
any angles of incidence [lOl [ll| . The incident field pumps 
those plasmon modes up and this lead to enhancement 
of evanescent waves reaching the surface on the image 
side of the slab further from the source. The induced 
displacement currents re-emit the light that reconstructs 
the image of the source without loss of resolution for 
features smaller than A. Similar effect involving surface 
plasmons is also responsible for extraordinary transmis- 
sion in thin metallic hole arrays [l^. This behavior is 
quite fragile, however, and limited by losses and spatial 
dispersion (metamaterials granularity) [l^ [l^, [l3| . Nev- 
ertheless, it has indeed been demonstrated experimen- 
tally by Lagarkov and Kissel that the NIM slab built 
with SRR interspersed with the wire mesh is able to 
resolve features ~ A/6 in the source separated by for 
microwave radiation (/ = 1.7GHz)[l^. Later on, the 
sub-wavelength resolution (~ A/6) was demonstrated by 
N. Fang et al. in the visible range using silver slab as 
a plasmonic medium ^l^. The silver slab is not a ma- 
terial with both permittivity and permeability negative 
{double negative), it only has e < as any metal at 
frequencies below plasmon frequency, and > 0, but 
the effect is still possible because the system is in qua- 
sistatic limit where the sign of permeability /i drops out 
of the result for image intensity [95. Various papers de- 
scribe metamaterials that show NIM-like response at far 
infrared frequencies p7| and, most interes ting ly, in near- 
optical and optical interval [ll [H, [13, [HI, [22[|. 



II. ELECTROMAGNETIC RESPONSE OF 
METALLIC NANOPARTICLE ASSEMBLIES 

The metallic periodic systems like 'fishnet' structure 
[Tl ] support an infinite set (bands) of electromagnetic 
(EM) waves, u = uj^, where k is the wavevector (quasi- 
momentum), n — 1,2,... the number of the band (Flo- 
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quet mode). At small wavevectors k such a crystal can 
be characterized by effective permittivity e and perme- 
ability n^^. Because of strong dispersion, it is easy to 
find crystals where some of the higher bands (second or 
third) support backward waves corresponding to nega- 
tive group velocity, see Refs.[3, S|- If such a band exists 
alone in a particular (usually narrow) frequency range, 
the crystal would operate as a NIM at those frequencies. 
Eleftheriades et al. have considered recently a possibility 
of an isotropic 3D NIM crystals (3D transmission lines) 
that might support backward waves in the first band if 
made of lumped capacitors and inductances, Ref. 24j, see 
Fig. la. This is directly follows from the 3D transmission 
line (TL) model suggested by Kron for Maxwell equations 
in isotropic space with positive permittivity and perme- 
ability, e,/Lt > 0^]. In the case of 'positive' medium the 
dispersion at small fc-vectors is positive for the lowest 
energy band, where eigenfrequency w — uJk is a growing 
function of fc. It should become a decreasing function of k 
after a dual transformation of the TL. Indeed, this is the 
standard passband-to-stopband transformation, where 
the dispersion remains exactly the same after a substi- 
tution C ^ L with obvious replacement uj ^ uj' ^ , 
which means that the group velocity for the dual band 
changes sign: Vg — du/dk > v'g{uj') < 0. The dual 
transformation does indeed result in negative dispersion 
in the doubly degenerate first band around the F point 
[2^ . However, the Kron TL lattice contains a few ele- 
ments per unit cell and, as a consequence, there is also an- 
other 'spurious' band present in the first Brillouin zone. 
The 3D TL crystal can be fairly well impedance matched 
to free space in ~GHz range and it supports a backward 
wave. The 'spurious' band can couple to highly atten- 
uated evanescent waves, but losses may be the limiting 
factor in subwavelength resolution experiment with dual 
Kron's TL lattice rather than the coupling to the 'spuri- 
ous' band [13]. 

It would be interesting to find an implementation of 
the dual Kron 3D transmission line model for NIM in 
visible range. With this in mind, Engheta et al. have 
speculated that Ag nanoparticles can play a role similar 
to lumped inductance at optical frequencies. One can 
try, for instance, substituting lumped inductances by Ag 
nanoparticles in the dual Kron lattice, hoping that the ca- 
pacitive coupling between the particles would bring about 
the same band structure as in 3D TL[2^, see Fig. lb. The 
nanoparticle system in Fig. lb is a poor representation 
of 3D TL system shown in Fig. la: although one does 
have both electric and magnetic response from metal- 
lic nanoparticles in ac field, they cannot be considered as 
lumped circuit elements. Indeed, the interaction between 
particles beyond nearest neighbors is likely important in 
open structures like Kron's cubic lattice, in addition to 
their non-negligible self-capacitance. We have calculated 
the two cases: (i) lattice of Ag spheres in vacuo (dielectric 
constant of the matrix is unity, rim = 1 , and (ii) dielec- 
tric constant of the matrix is = 1.4. The second case 
corresponds to experimentally accessible case of metallic 



spheres embedded in an organic matrix. We see from 
Fig. Id that the realistic case of Ag spheres in the dielec- 
tric matrix the electromagnetic response of the system is 
'soft' (Re(n) is small) but positive. Correspondingly, the 
dual Kron's lattice of Ag nanoparticles in vacuo shows 
only a weak negative index behavior, Re(n) w —0.3 but 
it is overwhelmed by losses, since Im(n) 3> |Re(n)|. It is 
likely, therefore, that the weak negative index behavior 
of Ag spheres in dual Kron's lattice won't be observable 
in sub-wavelength focusing experiments. 

To get more insight into electromagnetic (EM) re- 
sponse of nanoparticle arrays, we have performed 
extensive Finite Difference Time Domain (FDTD) 
modeling[23| of dual Kron lattices (Fig. 1) and close 
packed layers of spherical Ag nanoparticles (Figs. 2-4). 
The dielectric constant of Ag was assumed to have a 
Drude form 



with Up = 9.04 eV, 7 = 0.02 eV 28]. The bare Mie res- 
onance of Ag nanoparticles would be at ujm = ijJpl\f2> = 
5.22 eV or wavelength \m = 238 nmjlg. The trans- 
mission characteristics of the slabs have been calculated 
by FDTD and then used to estimate the effective index 
of refraction n, permittivity e and permeability /i from 
complex scattering coefficients according to a standard 
procedure ^30]. 

The calculated effective permittivity and permeabil- 
ity of the dual Kron's lattice of Ag particles with radius 
R — lOOnm show a series of sharp resonances, Re(e) < 
at A > O.S^m, and Re(/i) becomes negative at A < l/xm. 
In the same region Re(e) is also negative, and we see that 
the material has a small negative index at A w 0.9/im, 
where Re(n) ~ —0.3. However, the losses are very large 
there, Ini(n) w 1.7, and this will likely preclude sub- 
wavelength resolution with the lens made with this meta- 
material. Large losses seem to be quite general for sys- 
tems of metallic nanoparticles .26] and may severely limit 
their use for the purpose of sub- wavelength resolution. 

It is easily understood that the nanospheres should be 
almost touching to increase capacitive coupling between 
them. To see this, we considered Kron's lattice with Ag 
nanoparticles with radius R = 300nm with the gap of 
300nm between then, embedded in a material with dielec- 
tric constant em — 6.5. In this case the index remained 
positive at all wavelengths of interest, Re(n) « 2.5 at 
A = 1.2 — 3.2/im, since the displacement current were in- 
effective in producing Re(/i) < 0, and Re(e) only dipped 
below zero at a couple of very narrow intervals around 
resonant wavelengths. Since the volume fraction of the 
metal was quite small, the losses were minimal, typically 
Im(n) < 0.3. To get NIM behavior, one needs to sharply 
increase the metal fraction, but this immediately leads 
to prohibitively large losses, as described above. 

We have also studied arrays of Ag nanospheres with 
radius R = 30nm close packed into one, two, and three 
monolayers. Figs. 2-4. Such systems can be prepared, for 
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FIG. 1: Schematic of dual Kron lattice [24| (a) approxi- 
mated by the cubic lattice of Ag nanoparticles with radius 
7? = 100 nm and lattice period a — 300nm (b). Particles 
themselves show some inductance (marked by the coil sign) 
and are capacitively coupled to each other, as schematically 
indicated in the graph (b). There is some resemblance to 
the lattice of lumped elements[23] (a), (c) The effective per- 
mittivity e(aj) and permeability fJ,{u}) of a slab (2 u.c.) of Ag 
nanoparticles (b) . (d) Real and imaginary parts of index n, for 
Kron's lattice in vacuo and in matrix material with = 1-4. 
Ag lattice in the matrix shows some 'softness' at A « 0.9/im 
[small Re(n)], whereas system in vacuo exhibits a weak nega- 
tive index behavior with Re(n) « —0.3 but Im(n) 3> |Re(n)|. 



instance, by Langmuir-Blodgett technique. The particles 
have been stacked in a triangular close packed monolayer 
(A), close packed bilayer with AB packing (B is a trian- 
gular layer of nanospheres placed on top of the A layer so 
that the top sphere touches three underlying spheres in 
the layer A). Note that there are 'see-through' channels 
of deep sub- wavelength size in either A monolayer or AB 
bilayer that facilitate high transmission of at A 1.1/im, 
Figs. 2-4. The effect is apparently similar to extraordi- 
nary transmission of li ght through sub- wavelength hole 
arrays in metallic films [l^]. It is worth comparing this 
with Kron's system, Fig. lb, which is an open cubic lat- 
tice and has high transparency even when a slab with 3 
unit cells is considered. The transmission remains high 
even for touching or slightly overlapping nanoparticles, so 
the assembly is transparent and conductive at the same 
time. Quite expectedly, the transmission of AB bilayer is 
reduced in comparison with a monolayer, from T « 80% 
to about 60%, but it still remains very large. 

System with three layers of nanospheres packed in 
ABC sequence characteristic of face-centered cubic lat- 
tice (fee) is interesting, because there is no 'see-through' 
channels in ABC stack. However, even in this case trans- 
mission exceeds 40% at the resonance. Fig. 4. We see 
from Figs. 2-4 that one to three monolayers of close 
packed metallic Ag nanospheres produce an extraordi- 
nary transmission in the vicinity of A = 1.1/zm, which 
is much larger than the radius of the spheres R and the 
lattice spacing a. It is also much larger than the Mie res- 
onant wavelength Xm = 238 nm. Since in ABC trilayers 
there is no open "channels" for light to squeeze through 
the film, the explanation of unusual 'transparent metal' 
behavior lies in the fact that the incident light strongly 
couples to an array. Strong coupling of the incident light 
is obviously facilitated by periodic 'roughness' of the ar- 
rays, since the (quasi)momentum conservation is easier 
to meet. It would be interesting to see how surface plas- 
mon polaritons are supported by arrays of nanoparticles, 
in other words, what kind of surface plasmon waveguid- 
ing is possible with arrays of nanoparticles considered in 
the present paper [s^]. 

Importantly, electric field concentrates in the region 
where the spheres touch, especially when the centers of 
the spheres are oriented along a polarization of incident 
electric field, which looks analogous to corresponding 
electrostatic problem of two close metallic spheres. In 
fact, the local electric field enhancement exceeds a factor 
of 77 ~ 30, which would facilitate strong Raman signal if 
the species were positioned on the particles near the field 
peaks (Raman enhancement factor ry'* ~ 10^). 

In terms of electromagnetic response, there is a clear 
topological difference between a monolayer and a few- 
layer systems. Indeed, at normal incidence the magnetic 
field in the incident wave cannot excite any displacement 
currents that would for a close loops, since the field does 
not 'see' any such loops. In AB and ABC layers there 
are three-membered loops in the fee lattice that can pro- 
duce such a magnetic response. As follows from the data 
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FIG. 2: (Top panel) Transmission T, reflection R, and ab- 
sorption ^ of a monolayer of closely packed nanospheres with 
radius a = 30nm. Note a very large transmission T ~ 0.75 
at A=l.lum, much larger than the period of the monolayer, 
A ^ a. (Bottom panel) Permittivity e(a), permeability ^, and 
the effective index n of the monolayer of Ag nanoparticles. 



in Figs. 3,4, we do see some 'softness' in magnetic re- 
sponse of the AB and ABC arrays in Ke{fi) at A < 0.8/im, 
Fig. 2b. At the same time, the index for the monolayer is 
positive and does not show any particular 'softness', cf. 
Fig. 2. 



III. CONCLUSIONS 

We have analyzed various arrays of metallic nanopar- 
ticles in search for an isotropic negative index be- 
havior that was envisaged for dual Kron's cubic lat- 
tice of lumped circuit elements (3D transmission line 
model) [23] • In particular, we wanted to see if the assem- 
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FIG. 3: (Top panel) Transmission T, reflection R, and ab- 
sorption y4 of a two (AB packing sequence) layers of closely 
packed nanospheres with radius a = 30nm. Note a very large 
transmission at A=l.lum, much larger than the period of the 
monolayer, A ^ a. (Bottom panel) The effective index n of 
the bi-layer of Ag nanoparticles. We observe a transmission 
in excess of 50% at A=l.l/im. 



bly of metallic nanoparticles can exhibit effective nega- 
tive index behavior [2q. We have found that dual Kron's 
lattice of closely spaced Ag nanoparticles does not exhibit 
negative index behavior while embedded in a matrix, as 
we have shown for a matrix with rim = 1-4. Kron's lat- 
tice in vacuo may show a weak negative index behav- 
ior. However, it is overwhelmed by losses that would 
preclude sub-wavelength resolution with a slab of this 
metamaterial. We also made an extensive FDTD study 
of close-packed fcc-like arrays of Ag nanospheres. Ap- 
parently, they do support collective plasnion excitations 
and extraordinary transparency. Even in the opaque fcc- 
like trilayer of Ag nanospheres, where there are no see- 
through channels for light, the plasmons do transfer the 
excitation from the front to back surface and re-emit it, 
providing for high transparency in excess of 40-50%. In a 
monolayer of Ag particles the transparency is in excess of 
80%. Interestingly, the monolayer does not show electro- 
magnetic 'softness' in response, whereas in both bi- and 
tri- layers (AB and ABC packed films) we have found that 
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Re(n) « 0.3 at A « 0.7/xin. Although the arrays do not 
behave as negative index media, we see an intriguing be- 
havior characteristic of a 'transparent metal'. All those 
predictions would be very interesting to test experimen- 
tally. 




Wavelength (^m) 





FIG. 4: (Top panel) Transmission T, reflection R, and ab- 
sorption ^ of a three (ABC packing sequence) layers of closely 
packed nanospheres with radius a = 30nm. (Bottom panel) 
The effective index n of the bi-layer of Ag nanoparticles. Al- 
though there is no 'see-through' channels in this structure, we 
still observe a large transmission at \=\.\\xm. 
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